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Group-theoretical construction of extended baryon operators 
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The design and implementation of large sets of spatially extended baryon operators for use in lattice simulations 
are described. The operators are constructed to maximize overlaps with the low-lying states of interest, while 
minimizing the number of sources needed in computing the required quark propagators. 



1. OVERVIEW 

The Lattice Hadron Physics Collaboration 
(LHPC) was given a charge a few years ago by the 
late Nathan Isgur to extract as much as possible 
of the low-lying spectrum of baryon resonances 
from Monte Carlo simulations in lattice QCD. 
To accomplish this ambitious task, we have been 
pursuing two different approaches to construct- 
ing baryon interpolating field operators. Progress 
with the first of these approaches is reported else- 
where in this conference 021 • In this talk, de- 
velopments in the design and implementation of 
extended baryon operators using the second ap- 
proach are outlined. Initial efforts using this sec- 
ond method were described at the Tsukuba sym- 
posium last year [3]. Some details presented in the 
earlier reference will not be repeated in this brief 
report. 

The usual operator construction which mimics 
the approach one would take in continuous space- 
time (such as in the QCD sum-rule method) is 
very cumbersome, especially when tackling an en- 
tire spectrum, as opposed to a single state or a 
few states. Our approach described here is in- 



stead to directly combine the physical charac- 
teristics of baryons with the symmetries of the 
lattice regularization of QCD used in simula- 
tions. For baryons at rest, our operators are 
formed using group-theoretical projections onto 
the irreducible representations (irreps) of the Oh 
symmetry group of a three-dimensional cubic 
lattice. There are four two-dimensional irreps 
Gig, Giu, G2g, and two four-dimensional rep- 
resentations Hg and The continuum- limit 
spins J of our states must be deduced by exam- 
ining degeneracy patterns across the different Oh 
irreps |3j. 

Baryons are expected to be rather large ob- 
jects, and hence, local operators will not suffice. 
Our approach to constructing extended operators 
is to use covariant displacements of the quark 
fields. Displacements in different directions are 
used to build up the appropriate orbital struc- 
ture, and displacements of different lengths can 
build up the needed radial structure. All of our 
three-quark baryon operators are linear superpo- 
sitions of gauge-invariant terms of the form 



ABC 



where A, B, C indicate quark flavor, a, b, c are 
color indices, a, /?, 7 are Dirac spin indices, ip indi- 
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cates a smeared quark field, and denotes the 
p-link covariant displacement operator in the j-th 
direction. The three-dimensional gauge-covariant 
Laplacian is used to smear the quark fields: 

4'{x) = [l + gKr^ ^j{x), (2) 
A,p{x) = (t>fc(x)V'(x+fc)-V(a;)),(3) 

fc=±l,±2,±3 

where g and integer ng are tunable parameters, 
and Ufj_{x) denotes a link variable smeared us- 
ing the stout link algorithm 4 . The p-link gauge- 
covariant displacement operator is defined by 

D^Oix) = U,{x) . . . U,{x+(p-l)j)0{x+pj).{4) 

In the definition of ^^t^j-yk above, the dis- 
placement directions i,j,k can take values 
0, ±1,±2, ±3, where the value indicates that 
no displacement occurs. To simplify matters, the 
displacement length p is taken to be the same 
for all displaced quarks in a given term. Ul- 
timately, there are six different spatial orienta- 
tions that we use, shown in Tabled The singly- 
displaced operators are meant to mock up a 
diquark-quark coupling, and the doubly-displaced 
and triply-displaced operators are chosen since 
they favor the A-flux and F-flux configurations, 
respectively. 

Next, we combine the ^ai^pj-^k i^^^o elemental 
operators Bf {t, x) having the appropriate flavor 
structure characterized by isospin, strangeness, 
etc. We work in the to,i = nid (equal u and d 
quark masses) approximation, and thus, require 
that the elemental operators have deflnite isospin, 
that is, they satisfy appropriate commutation 
relations with the isospin operators T3,r-|_,r_. 
Since we plan to compute full correlation matri- 
ces, we need not be concerned with forming op- 
erators according to an SU (3) flavor symmetry. 
For zero momentum states, we impose translation 
invariance: Bf{t) = ^xBf{t,x). Maple code 
which manipulates Grassmann fields was used to 
identify maximal sets of linearly independent ele- 
mental operators. The numbers of such operators 
are listed in Table |21 

The final step in our operator construction is 
to apply group-theoretical projections to obtain 



Table 1 

Six types of three-quark ^^t^j-yk operators used. 
Smeared quarks fields are shown by solid circles, 
line segments indicate covariant displacements, 
and each hollow circle indicates the location of 
a color Eabc coupling. For simplicity, all displace- 
ments have the same length in an operator. 
Operator type Displacement indices 



i = j = k = 
single-site 

singly-displaced 

* ® * i = 7 = — fc fc ^ 

doubly-displaced-I ' ' 

« = 0, |j|v^|fc|, jk^O 
i = -j, \j\^\k\, jk^O 
\^\^\j\^\k\, ijk^O 



doubly-displaced-L 
triply-displaced-T 
triply-displaced-0 



operators which transform irreducibly under all 
lattice rotation and reflection symmetries: 

B^^it) = ^EneoA'^x {RrUj,Bf{t)ul (5) 
yOh 

where A refers to an Oh irrep, A is the irrep row, 
gOh is the number of elements in Oh, d\ is the 
dimension of the A irrep, Dmn{R) is a A repre- 
sentation matrix corresponding to group element 
R, and Ur is the quantum operator which imple- 
ments the symmetry operations. 

The projections in Eq. |SJ| are carried out using 
Maple in the following sequence of steps, (a) A set 
of Mb linearly-independent elemental operators 
B^ (t) that transform among one another under 
Oh is identified, (b) We obtain the Mb x Mb 
representation matrices Dij (R) which satisfy 

[/« Bf{t) Ul = EfJi Bf{t) Dy,{R). (6) 

(c) Since the resulting representations may not 
be unitary, we need to compute the Hermitian 
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Table 2 

Numbers of linearly independent elemental oper- 
ators. A^'^ operators have uuu flavor structure, 
(f2~ have sss), N'^ operators are uud — duu, S+ 
are uus and suu (EP are ssu and uss), and A*^ 
operators have uds — sud flavor structure. 



Operator Type 


A,n 


N 


E,S 


A 


single-site 


20 


20 


40 


24 


singly-displaced 


240 


384 


624 


528 


doubly-displaced-I 


192 


384 


576 


576 


doubly-displaced-L 


768 


1536 


2304 


2304 


triply-displaced-T 


768 


1536 


2304 


2304 


triply-displaced- 


512 


1024 


1536 


1536 



Table 3 

Total numbers of baryon operators for two differ- 
ent displacement lengths. 



Irrep 


\Vl 


N 




A 


Gig 


221 


443 


664 


656 


Giu 


221 


443 


664 


656 


G2g 


188 


376 


564 


556 


G2u 


188 


376 


564 


556 




418 


809 


1227 


1209 


Hu 


418 


809 


1227 


1209 



displacement orientations of the baryon sources. 
We hope to report on initial quenched simula- 
tions in the near future. Due to the large sizes 
of the baryons and their large masses, the use of 
improved actions on anisotropic lattices is impor- 
tant. Ultimately, when quark loops are included 
at realistically light quark masses, multi-hadron 
(baryon-meson) operators must be included in 
our correlation matrices, and finite-volume tech- 
niques will need to be employed to ferret out the 
baryon resonances from uninteresting scattering 
states. We are currently exploring different ways 
of building such operators. 

Although we focused on tri-quark baryon oper- 
ators, our approach to constructing hadronic op- 
erators is applicable to mesons, tetra-quark, and 
penta-quark operators. By including E and B 
fields into our operators, hadrons bound by an 
excited gluon field can also be studied. This work 
was supported by the U.S. National Science Foun- 
dation through grants PHY-0099450 and PHY- 
0300065, and by the U.S. Department of Energy 
under contracts DE-AC05-84ER40150 and DE- 
FG02-93ER-40762. 
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metric matrix M 

M., = — Efleo- Efe=i D,,{RrDk,{R). (7) 

(d) For each irrep A, we compute the Mb x Mb 
projection matrix for row A = 1: 



9oP " 



(e) From the rows of this projection matrix, r 
linearly-independent operators are obtained: 



B^'''{t)=E-=i 4'^ Bfit), (A = l) (9) 

Em=i 4'^* cf,'^ = (z = 1. . .r). (10) 

(f) For each of the r operators B^^^{t) in the 
first row A = 1, we obtain partner operators in 
all other rows fi > 1 using 



Mb , 



M 

A^F _ AAF 

We average over rows for increased statistics. The 
total numbers of operators one obtains in each 
row in the different isospin channels are given in 
Table El 

With such large numbers of operators, mini- 
mizing the number of quark propagator sources 
is crucial. We have recently finished writing soft- 
ware which, for each element of the baryon corre- 
lation matrices, performs all of the needed Wick 
contractions and simplifies the results by rotat- 
ing individual terms to minimize the number of 
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